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We used the atomic force microscope to manipulate and unfold individual molecules of the titin 
127 domain and reconstructed its free energy surface using Jarzynski's equality. The free energy 
surface for both stretching and unfolding was reconstructed using an exact formula that relates 
the nonequilibrium work fluctuations to the molecular free energy. In addition, the unfolding free 
energy barrier, i.e. the activation energy, was directly obtained from experimental data for the first 
time. This work demonstrates that Jarzynski's equality can be used to analyze nonequilibrium 
single-molecule experiments, and to obtain the free energy surfaces for molecular systems, including 
interactions for which only nonequilibrium work can be measured. 

PACS numbers: 87.15.He, 87.14.Ee, 87.64.Dz 
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One way to probe molecular properties is to drive a 
system out of equilibrium and to observe the response. 
Interpretation of data from dynamic measurements al- 
lows one to reconstruct both the equilibrium proper- 
ties of molecules and responses to external perturbations 
1] . Equilibrium parameters are usually deduced from ki- 
netic measurements, and it remains challenging to relate 
nonequilibrium distribution data to equilibrium proper- 
ties [2|. Advances in single- molecule manipulation and 
measurement techniques have made it possible to directly 
probe the dynamics of molecular interactions 0, @| . The 
noncquilibrium work theorem, i.e. Jarzynski's equality 
d, relates nonequilibrium measurements of nanoscale 
systems to equilibrium properties [1, 0, S] • It promises to 
extract thermodynamic parameters such as free energies 
from single- molecule measurements. 

Forced unfolding of single molecules, now achievable 
using the atomic force microscope (AFM) and laser opti- 
cal tweezers, has been used to probe the molecular inter- 
actions and mechanical properties of individual molecules 
d, [§]. In these experiments, single molecules are held 
at both ends and stretched while the cantilever spring 
restoring force (F s ) is measured. Applying an external 
force drives the system out of equilibrium, and transi- 
tions between states are directly observed as the system 
settles to a new equilibrium state. However, interpreta- 
tion of these results and deduction of equilibrium prop- 
erties from these nonequilibrium measurements remains 
controversial 0, El EE El • 

It has been widely anticipated that equilibrium free 
energy differences can be derived from nonequilibrium 
measurements using Jarzynski's equality 5]. The differ- 
ence in equilibrium free energy, AG, is related to the 
fluctuations of work_performed during a nonequilibrium 
process, W\, by 0, 

(e-^) N = I dWxp(W x )e-P w * = e^ AG (1) 



resents an average over N realizations of the process, and 
the equality is exact in the limit N — > oo. The nonequi- 
librium work distribution, p(W\), depends on the sched- 
ule for varying the work parameter A, which is the ex- 
ternal perturbation. The equality is simple; however, its 
application to interpreting single-molecule results is not 
straightforward. The equation involves the thermody- 
namic work done on the system and the controlled work 
parameter with W\ = J F-dX. In AFM experiments, the 
system includes the cantilever spring and the molecule 
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where /3 = (1/fcgT), ks is the Boltznrann constant, and 
T is the temperature of the thermal bath. The (...) jy rep- 
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FIG. 1: (color). Single-molecule pulling experiments using 
AFM. (a) One end of the molecule is attached to the cantilever 
tip and the other end to a gold substrate, whose position is 
controlled by a piezoelectric actuator. An analogue of the 
single-molecule force measurements is illustrated. The can- 
tilever spring obeys Hooke's law, whereas the protein molec- 
ular spring follows the worm-like chain model (illustrated us- 
ing rubber bands), (b) A representative force versus time 
trace, taken at 1.00 /im/s using a cantilever with a spring 
constant of 0.04 N/m. Each force peak represents unfolding 
of an individual titin 127 domain, with the final peak result- 
ing from the detachment of the molecule from the AFM tip. 
(c) Corresponding force-extension curve. The tip force base- 
line was determined using the part of the force curve where 
the molecule is completely detached from the tip, when the 
cantilever spring is at its equilibrium position. 



1 



PHYSICAL REVIEW LETTERS 



(2007) to appear 



plus water, and A refers to the change in cantilever an- 
chor to stage distance (see Fig. [1]) , not the tip-to-sample 
distance, which measures the molecule end-to-end dis- 
tance z, i.e. the order parameter, or reaction coordinate. 

We briefly review the experimental setup to which 
Eq. dl]) applies. Consider at time t = 0, the system is 
at an equilibrium state A(0) = Xa- We perform external 
work on the system by controlling the work parameter 
following a pre — determined schedule, A(i), from an ini- 
tial state Xa to a final state Xb- The system is then 
allowed to relax to equilibrium while A is held constant 
at As. Since we do not perform external work on the 
system during relaxation, we can omit this last step and 
the equality still holds. Hence Jarzynski's equality allows 
us to determine the G(A) from an equilibrated state A to 
an arbitrary state B. 

A proof-of-principle experiment and molecular dynam- 
ics simulations testing the Jarzynski estimator have been 
performed 
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1 1 51 ] . The experimental test involved 
stretching individual RNA molecules reversibly and ir- 
reversibly using optical tweezers, and the free energy of 
unfolding, i.e. the stability of the molecules, was deter- 
mined. However, the usefulness of Jarzynski's equality 
lies with its ability to obtain directly the entire free en- 
ergy landscape, which could only be estimated using ki- 
netic approaches to date [l6[ . We will show that Jarzyn- 
ski's equality can be used to determine directly the free 
energy profile of molecular stretching and unfolding, in- 
cluding the free energy barrier of unfolding. 

Our system of interest is the mechanical unfolding of 
the 127 domain of human cardiac titin [17]. The me- 
chanical properties of the immunoglobulin (Ig)-like do- 
mains are directly correlated with the protein's biolog- 
ical function in the muscles 0]. The kinetic barrier of 
these mechanical proteins is important in determining 
the dynamic behavior of the proteins during the stretch- 
release process. Therefore, the titin free energy surface, 
including the unfolding barrier height, is useful for quan- 
tification of titin's function in the heart muscle. 

We used AFM to stretch individual molecules of eight 
serially linked repeats of the titin 127 domain, as illus- 
trated in Fig. Q] The protein was stretched when the 
substrate stage was moved by A, which was set at a con- 
stant velocity v, i.e. X = vt. The cantilever displacement 
from its equilibrium position Az was recorded, and the 
molecular end-to-end distance as a function of time was 
calculated using z — A — Az. The force curves are aligned 
using the best worm-like chain (WLC) fit of the force be- 
low the unfolding force. This method has been shown 
to minimize the effect of instrument drift that affects 
the measured values [HJ . To correctly calculate AG as a 
function of the molecular end-to-end distance, we used an 
exact expression that connects the nonequilibrium fluc- 



tuations of work to the Gibbs free energy G(z) 0] 



where zq and z t are the end-to-end distances of the 
molecule at times and t during one realization of the 
process, F m is the force on the molecule, W z (t) is the me- 
chanical work done on the molecule up to time t, S(z — z t ) 
is the Dirac 5 function, and Uq is the potential energy 
stored in the cantilever spring at time 0. 

To calculate G(z) using Eq. ©, we divided each of the 
N trajectories of duration r into discrete time steps St 
so T = r/St, where T is the total number of time steps 
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FIG. 2: (color). Free energy reconstruction of titin 127 for 
pulling velocities of 0.05, 0.10, and 1.00 /^m/s obtained using 
64, 132, and 226 curves, respectively, (a) Typical unfold- 
ing force versus time curves for titin 127 domain taken at 
1.00 /im/s. Shown are 20 curves smoothed using a smoothing 
spline for display purposes, (b) Free energy G(z) calculated 
using the Jarzynski estimator, Gje, applied to the raw data. 
The averaged work, (W z ) = J^f W„/N, where W z = J F-dz, 
is displayed for comparison. (W z ) is larger than the equilib- 
rium free energy Gje by about a factor of 2 and is velocity 
dependent, whereas Gje is velocity independent. The curves 
are accurate up to the transition state (solid line), (c) Dis- 
tributions of work for z as a function of pulling velocity. The 
calculated work includes stretching and unfolding one domain. 
The curve fit to each distribution is a smoothing spline fit to 
the data as a guide to the eye. 
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in a given trajectory, 

ex P R?G(z(™))] « ^£$>( 



AT T 

Tl — 1 S=l 

exp(-/?[W n , s - U{z n , , A A )]X3) 



where z n>s is z at the s'th time step for the n'th trajec- 
tory, z n f) is the initial value of z for the n'th trajectory, 
and W n ,s is the work performed up to time t s — s ■ St for 
the n'th trajectory. We divided the z-axis into bins of 
width e and let z^ m > represent the mid-point of the m'th 
bin. The S function is l/e when z n>t falls inside the m'th 
bin and otherwise. The integration starts from the be- 
ginning of the curve, where the cantilever is close to its 
resting position, z = 0, at t = 0. This initial condition is 
required for using Jarzynski estimator, which states that 
the process needs to start from an equilibrated state. It 
is also advantageous when using Eq. ([3]) that the initial 
energy stored in the cantilever spring, Uq(z q , Xa), is close 
to 0. We compare this result to the approximate free en- 
ergy surface derived from 
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The unfolding free energy surface of titin 127 determined 
from Eq. (J3J) and Eq. © are very similar, perhaps due 
to the relatively stiff cantilever used in AFM. However, 
it is physically and theoretically more meaningful to use 
Eq. ^ , since determination of the entire free energy sur- 
face relies on converting the coordinate from t to z. Fig- 
ure [2] displays the free energy surface measured at three 
different velocities, determined using Eq. ([3]). 

The unfolding free energy barrier AG* can be calcu- 
lated from the reconstructed free energy curve. Using 
0.6 nm as the distance between the native and the tran- 
sition state (x u ) 11, HI, we calculated the unfolding free 
energy barrier AG* for pulling velocities of 0.05, 0.10, 
and 1.00 (im/s, to be 11.0, 11.7, and 11.4 kcal/mol, from 
Eq. © and 11.5, 11.5, and 10.7 kcal/mol from Eq. Q, 
respectively. The uncertainty in the averaged AG* = 
11.4 and 11.2 kcal/mol, calculated using the bootstrap 
method, is 0.4 and 0.3 kcal/mol for Eqs. © and g|, 
respectively. This result compares favorably to an es- 
timated value of 10-16 kcal/mol [IH, [2^, 21 1. A major 



source of error for AG* from Jarzynski estimator comes 
from the uncertainty in x u . Using the largest estimated 
error of 0.07 nm uncertainty in x u [HI], the estimated 
uncertainty of AG* is 1.2 kcal/mol. 

The free energy surface is accurately reconstructed 
from z = to 17 nm, the transition state. The free energy 
of unfolding is insensitive to the distance of reconstruc- 
tion. As an example, if we use 15 nm or 19 nm, AG* 
changes by 2.5 and 0.4%, respectively, for the pulling ve- 
locity of 0.05 /im/s. The vast majority of the proteins 
in the ensemble are in the folded state (99.9997 % using 
the free energy AG U from Ref. 22j) so the contribution 



from the initially unfolded proteins is negligible. To mini- 
mize the contribution from other unfolded domains to the 
measured free energy, we analyzed only the first domain 
stretching event. Using all domain unfolding events in 
the analysis changes AG* by less than 2%. 

Note that it is not possible to compare our results 
directly to published values, since AG* has not been 
determined, and only kinetic information is available. 
Chemical denaturant studies gave an estimated unfold- 
ing rate constant, of 6xl0~ 4 s~ 1 21], while forced- 
unfolding studies gave an estimated fe° of 1Q -3 — 10 — 6 s — 1 

was then calculated using 
k® = fcoe~' 3AG ". Since the prefactor ko of protein un- 
folding is unknown, the free e nerg y barrier can only be 
estimated by this procedure [24, 25]. However, combining 
our free energy determination with the kinetic informa- 
tion, we can determine the prefacto r l/fcp to be 6 /is, 
which lies within the expected range [2J, [26( . 

The free energy surface immediately past the transi- 
tion state cannot be reconstructed with high accuracy 
from constant velocity unfolding experimental data. This 
is because the force exerted on the molecule is discon- 
tinuous when the domain ruptures and expands against 
the cantilever. In the region where the domain ruptures 
and the cantilever snaps back to its equilibrium posi- 
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FIG. 3: Free energy surface of titin 127. (a) A typical force 
versus extension curve. The gray curve is the WLC fit to the 
following domain. The shaded region indicates that, when 
the domain ruptures and the cantilever snaps, the force on 
the molecule is not registered and, therefore, the free energy 
surface may not be recovered with high accuracy. The dashed 
line indicates an approximation of the force exerted on the 
molecule, (b) The free energy surface of unfolding titin 127. 
The cantilever position and the molecular extension at each 
stage are illustrated. The curve is composed of the recon- 
structed free energy surface up to the transition state (solid) 
and estimated free energy change [2^] and distance [23|] be- 
yond the transition state (dotted). 
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tion, the assumption that the force on the molecule (F m ) 
is balanced by the cantilever spring restoring force (F s ) 
no longer holds. Therefore, using the measured F s gave 
rise to an overestimate of the free energy. Note that 
even though the snapping process is almost instantaneous 
(small change in t, hence A), the change in z is significant 
because Az — F s /k s , where k s is the cantilever spring 
constant (see Fig. [3]). A lower pulling velocity and larger 
spring constant will reduce the size of the snapping re- 
gion. However, we can estimate the folding free energy 
barrier, AG fj, from the equilibrium unfolding free energy 
determined from chemical denaturant studies Q . Using 
AG U = 7.5 kcal/mol, we obtained AG^ = 3.9 kcal/mol, 
in the expected range for titin 127. Figure [3] summarizes 
the reconstructed free energy surface and its relation to 
pulling experiments. 

One requirement for using Jarzynski's equality is that 
the schedule for varying the work parameter A must 
be pre-determined [H, which means that constant force 
ramp is not an appropriate schedule. A constant dF s /dt 
requires force feedback and, since the measured force F s 
fluctuates from one pull to another, the result is a dif- 
ferent schedule of A for each realization. On the other 
hand, the dynamic force spectroscopy method commonly 
used in AFM pulling of proteins is particularly suitable 
for such analysis because the schedule for pulling is pre- 
determined and remains the same for all experiments per- 
formed at the same velocity. 

Using nonequilibrium single-molecule measurements 
and Jarzynski's equality, we have reconstructed the free 
energy surface of both mechanical stretching and unfold- 
ing of the 127 domain of human cardiac titin. Since the 
profile is an equilibrium property, the reverse of the free 
energy surface of stretching is equivalent to that of pro- 
tein folding from an extended state. The unfolding free 
energy barrier and the prefactor were determined directly 
from experimental measurements without having to as- 
sume either a two- or a three- state model, which are 
major sources of error in the event of populated interme- 
diate states. In fact, with adequate resolution and accu- 
racy, an intermediate state should be directly resolvable 
in the free energy curve. The topography and the rough- 
ness of the folding free energy landscape can also be de- 
termined. Reconstruction of free energy surfaces directly 
from experimental data is valuable to obtain fundamental 
thermodynamic properties such as the free energy barrier 
of unfolding, to understand the mechanical properties of 
the molecule, and to compare with theory and simula- 



tion results [27J . With a complete characterization of the 
free energy surface of molecular processes, questions such 
as whether thermal, chemical, and mechanical unfolding 
probe the same process may be resolved. Moreover, since 
the free energy surface is determined in a particular envi- 
ronment, how the free energy surface changes with envi- 
ronmental parameters such as temperature, solution ionic 
concentration, and acidity may now be evaluated. Quan- 



tification of the molecular response to external parame- 
ters should lead to a better understanding of molecular 
behavior in the complex cellular environment. 
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